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CD ' ABSTRACT 

> 

Q\ i Context. The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 
QO ■ it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging 
ly-j Cherenkov telescope (MAGIC-I). 

. Aims. We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. 
Methods. We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and 
optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelenght data within the syn- 
chrotron self-compton (SSC) framework. 
1 , Results. The source showed intense and prolonged y-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the 
. level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XKT) 

■ and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral 
• i-H : energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. 

\ Conclusions. The main outcome of the study is that within the homogeneous model high Doppler factors (40 < 6 < 80) are needed to reproduce 
5^ , the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models 

■ could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour 
scale variability. 

Key words. Gamma-rays: observations — Radiation mechanisms: non-thermal — Galaxies: BL Lacertae objects — individual:Mrk 421 
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1. Introduction 

Blazars, a common term used for flat-spectrum radio quasars 
(FSRQ) and BL Lacertae objects, constitute the subclass of ac- 
tive galactic nuclei (AGN) that is most commonly detected in 
the Very High Energy (VHE, E > 100 GeV) y-ray band. In 
these sources the dominant radiation component originates in a 
relativistic jet pointed nearly toward the observer. The double- 
peaked spectral energy distribution (SED) of blazars is attributed 



to a population of relativistic electrons spiraling in the magnetic 
field of the jet. The low-energy peak is commonly thought to be 
causedby synchrotron emission, because of its spectrum and po- 
larization. The second, high energy peak is attributed to inverse 
Compton scatte ring of low-energy photons in leptonic accel- 
eration models (jMaraschi et all Il992t iDermer & Schlickeiserl 
11993b IB! oom & Marscher, 1996). Alternative models invoking 
a relevant contribution from accelerated hadrons can also suffi- 
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ciently describe the obse rved SE Ds and light curve s (iMannheiml 
ll993tlMiicke et alJl2003l but see lSikora et aD2009l on FSRQs). 

Blazars are highly variable in all wavebands and the relation 
between variability in different bands is a key element in distin- 
guishing between different models. For instance, homogeneous 
leptonic models predict correlated variability between, e.g., X 
rays and y rays, which is alread y observed i n high-frequency 
peaked BL Lacs (HBL; see e.g. lFossati et alj|2008l on Mrk 421 
itself). On the other hand, phenom ena such as the "orphan" 
flare from 1ES 1959+650 reported in iKrawczvnski et all (|2004|) 
or the ultrafast (~ hundreds of seconds) events occasionally ob- 
served in some sources (e.g. Ahar onian et all 120071: 'Albert et all 
2007b) are harder to explain whithin this frame. 

Among blazars, HBLs are the most often observed subsam- 
ple in the VHE domain, because the high-energy bump peaks at 
GeV-TeV energies, while the low-en ergy peak is located at UV 
to X-rays energies dPadovanil 12007b . This makes HBLs, such 
as Mrk 421, ideal targets for sensitive, low-energy threshold 
imaging air Cherenkov telescopes (IACT) such as MAGIC-I, 
in combination with soft X-ray telescopes, which observe the 
synchrotron bump instead; this combination of instruments 
samples the source SED and unravel the regions of the two 
peaks, the most valuable tracers of the source state. 



Mrk 4 21 is one of the closest (z = 0.031. Ide Vaucoleurs et al.l 
1 199 lh and brightest ext ragalactic TeV sour ces; therefore it was 
the first to be detected (iPunch et all [T992) and remains one of 
the best studied. The VHE integral flux can vary from a few 
tenths to a few Cr ab Units (e.g ., see iDonnarumma et al . 2009; 
iHsu et alJ j2009 or Pichel 2009), on time scales as short as 15 
minutes dGaidos et all 119961) . The vF(v) distribution of the 
emitted photons follows the standard "double-bumped" shape, 
but varies significantly from low-activity states to the most 
intense flares, on ti me scales that in X-rays can be of few hours 
(lUshio et alll2009l) . The low-energy bu mp peaks in the 0.1-10 
keV range (see e.g. iFossati et all 12008*1) . as is usual for HBLs; 
the maximum of the high-energy bump is usually found below 
100 GeV, but can also move around according to the state of 
the source, usual l y following a "harder-when-brighter" be havior 
(iKrennrich eTail l2002t IFossati et all 120081: lAcciari et alll201 it 
Albert et al., 2007a j) analogous to th a t traced by the X-ray 
emiss ion (e.g. iBrinkmann etakl 120051: iTramacere et all l2007l 
120091) . 

This peculiar SED shape favors multi-wavelength (MWL) 
studies that exploit the MAGIC-I sensitivity and low-energy 
threshold in VHE, and soft X-ray telescopes. MAGIC-I can de- 
tect Mrk 421 at the 5cr level with exposures as short as a few 
minutes, depending on the source brightness. In the soft X-ray 
domain the All-Sky Monitor (ASM) onboard the Rossi X-ray 
Timing Explorer (RXTE) can continuously provide a daily aver- 
aged flux, while the X-Ray Telescope (XRT) onboard the Swift 
satellite can observe the source with far better precision and en- 
ergy resolution in ~ 1 ks targeted exposures. From the obser- 
vation of the source spectrum in both X-ray and VHE a unique 
set of physical parameters that describe the source can be derived 
within a single-z one synchrotron self-C ompton (hereafter, SSC) 
model, following iTavecchio et al.l (Il998l) . An SSC modeling of 
the SED of Mrk 421 has been already performed in the past 



(see e.g. iBednarek & Protheroel Il997t ITavecchio et al.l 1* 1998; 



iMaraschi et al.lll999l or the more recent IFossati et al 1120081)7 

L ately automated y 2 minim ization procedures (iFinke et all 

2008; Mankuzhivi l et all I20TIT) have been applied. 

The main limitations to previous works came from the use 
of the former IACTs such as Whipple (e.g. IFossati et all [2008 ) 



or HEGRA (e.g. [Takahashi et al., 2000), which were character- 
ized by a higher energy threshold and worse sensitivity at VHE. 
This in turn led to poor sampling of the IC peak region, basically 
limited to the less informative, steeply decaying hard energy tail 
of the bump; moreover, integration over different nights of ob- 
servation was commonly needed to obtain a significant VHE 
spectrum, thus yayeraging out the SED evolution, such as in 
IMaraschi et all d 19921) . In 2008 Mrk 421 went through a long 
and intense outburst phase, characterized by VHE fluxes quite 
constantly above the Crab level and superimposed shorter and 
brighter flares; a remarkably dense follow-up of this evolution 
was possible in optical, X-rays and VHE with MAGIC-I. The 
main outcome of this campaign is that eight tightly contempo- 
rary observations of Mrk 42 1 in optical, X-rays and VHE y-rays 
of active states could be achieved. This allowed the reconstruc- 
tion and modeling of the optical-UV/X-ray/TeV MWL SED of 
Mrk 421 on short time scales (~ 1 hour), the main improvement 
in this work with respect to the past literature on the subject. A 
similar approach was fol lowed, for the same Mrk 421, by the 
VERITAS Collaboration ( lAcciari et all 1201 ll) . 

The paper is structured as follows: in Section [2] we report 
on the observations and data analysis; in Section [3] we report 
the VHE light curve and spectra, and complementary results in 
X rays and optical-UV band; in Section [4] we build the SED 
of Mrk 42 1 in the eight states for which a set of simultaneous 
MWL data was available, which we model in the framework of 
a standard one-zone SSC model; finally we discuss the results 
in Section[5] 

2. Observations and data analysis 

The timespan of the MAGIC VHE observations of Mrk 421 re- 
ported here began in December 2007 and ended in June 2008. 
Contemporary data from other instruments are considered for 
the multi-wavelength analysis, namely soft X-ray data from 
RXTE/ASM, and Swift/XRT, optical-UV data from Swift/UVOT 
and optical R-band data from the Tuorla Observatory. A sum- 
mary description of the instruments, the datasets and the analysis 
follows. 

2.1. MAGIC-I VHE observations 

MAGICO-I (formerly MAGIC) is an IACT located on the west- 
ern Canarian island of La Palma, at the Observatory of Roque de 
Los Muchachos (28.75° N, 17.89° W, 2225 m a.s.l.). With its tes- 
sellated parabolic mirror (D - 17 m, f/D = 1), it has been the 
largest single-dish IACT in operation from late 2004 until the 
advent in 2009 of MAGIC-II, a twin (but substantially improved 
in many respects) telescope; since then the two telescopes are 
operated as a stereo IACT system (MAGIC Stereo). Its 234 m 2 
surface allowed for the lowest energy threshold among IACT 
systems at that time: the trigger threshold of the telescope at 
the epoch of this campaign reached as low as 60 GeV for ob- 
servations close to the zenith in optimal conditions. A detailed 
description of the tele scope and data analysis can be found in 
dedicated papers (e .g. Baixera s"et al.ll2004t ICortina et al 1 12001 
Albe rt et al.l*200 8a). All MAGIC-I observations considered for 
the prese nt study were perfor med following a major hardware 
upgrade dGoebel et al.. 2008) that was completed in February 
2007, which enhanced the time sampling capability of the data 
acquisition (DAQ) from 300 MHz to 2 GHz. This allowed a bet- 
ter rejection of the night sky background (NSB) and introducing 
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new refined analysis techniques (Tesc aro et all 12008) based on 
the time properties of Cherenkov signals, which lowered the in- 
tegral sensitivity for point sources down to 1.6 % of the Crab 
Nebula fl ux (for a 5cr sign ificant detection in 50 hours, above 
280 GeV: lAiiu et al.Ll2009h . 

During the observation period considered here, MAGIC ob- 
served Mrk 421 for a total of 81 nights, with exposure times 
ranging from ~ 20 to ~ 240 minutes. This comprised both 
short untriggered observations, aimed to an unbiased sam pling 
of the source state (studied in detail in Wa gner et al.ll201 lb . and 
deeper extended observations of peculiar states, triggered either 
by the former or by external alerts from other bands. All obser- 
vations were perfo rmed in the false-source tracking ("wobble", 
iFomin et al.lll994|) mode. The method consists of alternatively 
tracking two positions in the sky that are symmetrical with re- 
spect to the source nominal position and 0.4° away from it. 

The MAGIC-I data were analyze d using the sta ndard anal- 
ysis c hain described in lAlbert et all (l2008allcT) and lAliu et all 
(2009). Preliminary quality checks were performed to exclude 
poor quality data, such as those owing to bad weather or oc- 
casional technical problems. Furthermore, the dataset was re- 
stricted to observations performed under dark conditions, and 
in the range of zenith angle ranging from ~ 5° at culmina- 
tion to 46°. A cleaning algorithm involving the time structure 
of the shower images was then applied, which further selected 
the events and removed the NSB contribution to the images. 
Surviving images were p arametrized in terms of the extended 
set of Hillas parameters dHillasL 1 1985b described in the men- 
tioned literature. To suppress the unwanted background show- 
ers produced by charged cosmic rays, a m ultivariate class ifica- 
tion method known as ran dom forest (RF, [B reiman 2001) was 
implemented and applied dAlbert et al.L l2008bh . An analogous 
procedure allowed the estimation of the energy of the primary 
y-rays. The signal extraction was performed by applying cuts 
in the Size, Hadronness and Alpha parameters described in 
the aforementioned literature. In particular the Size cut, which 
we set to reject events with less than 150 photoelectrons of total 
charge, implied an energy threshold ~ 140 GeV in the present 
analysis. A total excess of ~ 48 x 10 3 events from the selected 
~ 60 hours of observation was detected. The whole analysis pro- 
cedure was validated step by step on compatible datasets from 
observations of the Crab Nebula. 

2.2. Optical, UV and X-ray observations 

The Swift satellite dGehrels et all |2004|) is a NASA mission, 
launched in 2003, devoted to observations of fast transients, 
namely prompt and afterglow emission of gamma-ray bursts. 
These are detecte d with the monitoring coded mask Burst Alert 
Telescope (BAT, Barthe Tmv et al.l 120051) which is sensitive to 
15-150 keV X-rays and covers a wide field of view (FoV) 
with a resolution of few arcminutes, and then rapidly targeted 
with t he two co-align ed pointing instruments, X-Ray Telescope 
(XRT. lBurrows et al.ll2005l) and Ultra- Violet Optical Telescope 
(UVOT lRoming et aUl2005h . 

The fast repositioning capability of the spacecraft allows 
snapshots of variable sources with little overheads. For Mrk 421, 
observations lasting ~ 1 ks allow the derivation of a detailed X- 
ray spectrum and multi-filter optical-UV photometry because of 
the sensitivity of the targeted instruments and the brightness of 
the source. 

Swift/XRT is a Wolter type-I grazing incidence telescope, 
with 110 cm 2 effective area, 23.6' FoV and 15" angular res- 
olution, sensitive in the 0.2-10 keV energy band. During the 



MAGIC campaign the instrument performed 43 targeted X- 
ray observations of Mrk 421 of typical exposure times 1— 

2 ks. Swift- XRT data were reduced using the software dis- 
tributed with the heasoft 6.3.2 package by the NASA High 
Energy Astrophysics Archive Research Center (HEASARC). 
The xrtpipeline was set for the photon counting or window 
timing modes and single pixel events (grade 0) were selected. 

UVOT is a 30 cm diffraction-limited optical-UV telescope, 
equipped with six different filters, sensitive in the 1700-6500 
wavelength range, in a 17' x 17' FoV. Unfortunately, during the 
January 2008 campaign UVOT did not observe the source, so 
that the UVOT datasets were fewer than the XRT pointings and 
therefire five out of the eight datasets studied in Section |4]i have 
no contemporary UVOT observations. Therefore we restricted 
the analysis of UVOT data to the three observations simultane- 
ous with MAGIC-I performed on February 1 1 and April 2 and 

3 with the UV filters alone. The analysis was performed with 
the uvotimsum and uvotsource tasks with a source region 
of 5", while the background was extracted from a source-free 
circular region with radius equal to 50" (it was not possible to 
use an annular region because of a nearby source). The extracted 
magni tudes were corrected f or Galactic extinction using the val- 
ues oflSchlegel et"ai1 (1 19981) and applying the formulae by iPeil 
(1992]) for the UV filters, and eve ntually were converted into 
fluxes following lPoole et alj d2008h . 

The All-Sky M onitor (ASM) on board the Rossi X-ray Timing 
Explorer (RXTE, B radt et al.l 1 993!) is sensitive enough to set one 
point per day from Mrk 421, which means a poorer precision, but 
denser coverage than Swift/XRT; 

The publicly available ASM data products were taken from 
the results provided by the ASM/RXTE teams at MIT and at the 
RXTE SOF and GOF at NASA's GSFC. 

The Tuorla Observatory constantly monitors the MAGIC 
VHE (known or potential) target sources, with the 35 cm re- 
motely operated Kungliga Vetenskaplika Academy (KVA) op- 
tical telescope that is also located at Roque de los Muchachos 
and with a 103 cm telescope located at Tuorla, Finland. During 
the period of the MAGIC-I observations, 117 photometric mea- 
surements of Mrk 421 were obtained in the Johnson 7?-band. 
The optic al data were reduced by the Tuorla Observatory as de- 
scribed in Nils son et al.l (120071) . The light contribution from the 
host galaxy and nearby companion galaxy (Fi, +Cg = 8.07 + 0.47 
mJy) was subtracted from the measured fluxes. 



3. Results 

3. 1. MAGIC-I VHE light curves 

The night-averaged integral flux above a conservative thresh- 
old of 200 GeV was calculated for each of the 66 nights with 
datasets that survived the quality cuts. The VHE light curve 
of Mrk 421 along the campaign is plotted in the top panel of 
Figure Q] Interestingly, al though Mrk 421 is be lieved to emit 
a low VHE flux baseline ( Schub nell et all 1 1996b . the flux was 
seldom below one Crab Unit (hereafter C.U. corresponding to 
an integral flux Fg> 2ooGev = 1-96 ± 0.05 stat X 10 cnT 2 s _1 , 
Albert et al. 2008a) for the entire period, confirming an in- 
tense and persistent active state. The maximum observed flux 
(f£>200Gev = 6.99 + 0.15 JM/ x lO-^cnrV 1 ) was on 2008 
March 30. Similar fluxes were reached in the flare that oc- 
curred in June 2008, whi ch was already studied in detail in 
iDonnarumma et all d2009b . Similar high flux le vels during the 
same period are reported in lAcciari et al.l (1201 lb . with a record 
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Fig. 1. Multi-wavelength light curves of Mrk 421 along the MAGIC-I observation period; full circles mark the fluxes observed 
when MAGIC-I and Swift/XRT were pointed at the source simultaneously. Upper panel: MAGIC-I VHE light curve above 200 
GeV, for the 66 observation nights that passed quality cuts. MAGIC-I detected the source clearly in all nights; the integral flux 
was below the Crab Unit (C.U., * f£>200GeV = 2.0 x 10~ 10 photons •cmT 1 s^ 1 , represented here by the dashed horizontal line) 
only in a few nights. A maximum flux of ~ 3.6 C.U. was observed on 2008 March the 30th (MJD=54555). Middle-upper panel: 
soft X-ray (0.2-10 keV) count rates measured by Swift/XRT. Middle-lower panel: soft X-ray (2-10 keV) count rates measured by 
RXTE/ASM. Lower panel: Johnson R-band optical light curve from the Tuorla Observatory. 



flux of ~ 12 x 10 "cm 2 s 1 above 300 GeV (corresponding to 
~ 10 C.U.) observed in May 2008. 

The sensitivity of MAGIC-I allowed us to investigate the 
sub-hour scale evolution of VHE flux for Mrk 421 in high state, 
searching for the ra pid variations already reported in literature 
(Gai dos et all fl996). The most interesting result was obtained 
on 2008 February 6, when a long (~ 4 hours) observation of a 



high (~ 2.5 C.U. above 200 GeV) state was performed. The VHE 
light curve in 8-minute time bins is shown in Figure [2] above a 
softer (E > 200 GeV, upper panel) and harder (E > 400 GeV, 
lower panel) energy threshold. An episode of variability with 
doubling/halving times down to 16 minutes can be seen with 
the harder cut. The hypothesis of a steady flux is unfit in both 
light curves according to results of a x 2 test, giving ^ 2 /d.o.f. of 
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55/28 (probability below 0.2%) and 63/28 (probability ~ 10" 4 ), 
respectively, which confirms variability on the scale of hours or 
less. Unfortunately, no simultaneous Swift/XRT observation was 
performed in this night, therefore we could not add this VHE 
observation to the set of simultaneous MWL SED. No firm con- 
clusion could be drawn on the sub-hour variability for the si- 
multaneous datasets, because some observation windows were 
very short (e.g. April the 2nd and 3rd) and other were made in 
lower flux levels (e.g. January 8, 9, 10), which led to poorer event 
statistics. 

3.1.1. Multi-wavelength data 

The Swift-XRT observed count rates in the 0.2 - 10 keV band 
are reported in the middle-upper panel of Figure [T] 

The count rates observed by RXTE/ASM in the 2 - 10 keV 
band are shown in the middle-lower panel of Figure Q] 

The R-band optical light curve from KVA observations is re- 
ported in the bottom panel of Figure Q] while the available mea- 
surements related to the simultaneous datasets listed in Table 
[3] (see Section HJ are plotted in Figure [6] after correction for 
Galactic extinction, again applied according to the values of 
ISchlegeletal] (119981) . 



MAGIC-I. Anyway, the VHE spectrum derived from this obser- 
vation is intriguingly hard, peaking around 500 GeV, well within 
the MAGIC-I band. In Figure [3] the observed SED (black open 
triangles) is plotted together with the deabsorbed one (red filled 
circles), which peaks above 1 TeV. 
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3.2. VHE spectra derived from MAGIC-I data 

We restricted the study of the spectra to the subset of the 
eight observations of interest for modeling the MWL SED (see 
Section |4]i, listed in Table [3] From each observation we derived 
a VHE spectrum in bins of the estimated energy of the y-ray 
prima ry events. Then we a pplied the Tikhonov unfolding algo- 
rithm (Al bert et all |2007cT) to reconstruct the physical spectrum 
in terms of the true energy of the primary y-rays. A best fit to 
the data was then performed, assuming a log-parabolic model 
for the differential spectrum: 



dN E .( fl+ £.iog(J-)) 

dEdAdt J y E 



(1) 



where the pivot energy Eq is chosen 300 GeV in the present case. 
In three cases a simple power law was sufficient to fit the data. 

For each night, the integral VHE flux above 200 GeV, the 
parameters of the fit to the observed (no EBL correction) emis- 
sion and the x 2 are reported in Table Q] quoted uncertainties are 
statistical only. 

Adopting the a parameter, giving the slope of the spectrum 
at the pivot energy, as an estimator of its hardness, it is evi- 
dent that the well-observed "harder when b righter" trend (see 
e.g. lAcciari et all 1201 it iFossati etafl 120081) is nicely repro- 
duced in these spectra. The spectral points were subsequently 
corre cted for extra-balactic ba ckground light (EBL) absorption. 
The lFranceschini et all (2008) EBL model has been assumed in 
this work, even if consistent re sults can be obtained w ith other 
models such as the more recent Domm guez et all (1201 lb model, 
given that for this close-by source the model-to-model differ- 
ences in opacity below 10 TeV are dominated by the statisti- 
cal uncertainties in hour-scale integrated VHE spectra. The data 
points are plotted in Figure [4] along with the SSC models (see 
Section |4~TI ) that are anticipated here only as a help to guide the 
eye. For comparison we also plot, without a model, the SED 
built from the observation achieving the highest VHE flux of 
the whole campaign (March 30). Unfortunately, we were un- 
able to include this interesting dataset in the SED study, because 
Swift could only observe with 14 hours of delay with respect to 



Fig. 3. VHE SED of Mrk 421 derived from the MAGIC-I ob- 
servations performed on 2008 March 30, when the flux rose to 
3.6 Crab Units. Data points before (black open triangles) and af- 
ter (red fi lled circles) applying a c orrection for EBL absorption 
following iFranceschini et al . 2008 are shown. The observed po- 
sition of the IC peak is evaluated at ~ 500 GeV from the fit with 
a curved power law, and above 1 TeV after deabsorption. This 
VHE spectrum was the hardest among the ones studied here, as 
illustrated in Figure |4] 



3.3. Soft X-ray spectra derived from Swift/XRT data 

For the eight simultaneous observations with MAGIC-I listed in 
Table [3] we extracted the Swift/XRT spectra to build the MWL 
SED (see Section 0]). Data were rebinned to obtain at least 30 
counts per energy bin. Broken power-law models were fitted to 
the spectra in the range 0.35-10 keV. The X-ray reddening due 
to absorbing systems along the light travel path was corrected 
assuming the Galactic value for th e column density of neutra l 
hydrogen N H = 1.6 x 10 20 ctrT 2 (lLockman & Savagel 119951) . 
Table [2] reports for each dataset (with uncertainties in parenthe- 
ses) the Obs ID, the UTC time at the beginning of observation, 
the exposure time, the integral flux in the 2-10 keV band, the 
spectral indexes, break energy and normalization at 1 keV of the 
broken power law, and the resulting x 1 (reduced x 1 ) with the 
number of degrees of freedom. 



4. Simultaneous multi-wavelength datasets 

Below we focus on the eight cases for which tightly simultane- 
ous observations in VHE with MAGIC-I telescope and in X- 
rays with Swift/XRT could be performed. Table [3] summarizes 
the observation logs of the two instruments for these nights. For 
each one we report the beginning and the end of the MAGIC- 
I observation time span, the total effective time, and the ZA 
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Fig. 2. Mrk 421 VHE light curves in 8-minute time bins from the observations taken on 2008 February 6. Integral flux of excess 
(filled circles) and background (thin crosses) events are plotted. The energy threshold is 200 Ge V (upper panel) and 400 GeV (lower 
panel). 



Night 


Integral flux 


/o 


a 


b 


X 2 /d.o.f. 




-2 -1 

cm s 


cm 4 s _l TeV~' 








yyyy-mm-dd 


[xlO^ 10 ] 


[xlO- 10 ] 










(E > 200 GeV) 


(E = 300 GeV) 








2008-01-08 


2.13 ±0.20 


5.9 ±0.7 


-2.72 ±0.12 




0.40/5 


2008-01-09 


2.61 ±0.11 


6.3 ±0.3 


-2.50 ± 0.07 


-0.44 ±0.15 


0.63/7 


2008-01-10 


2.53 ±0.16 


7.4 ±0.5 


-2.42 ± 0.08 


-0.52 ± 0.20 


1.26/6 


2008-01-16 


4.42 ±0.14 


10 ±1 


-2.25 ± 0.07 


-0.33 ±0.10 


0.34/6 


2008-01-17 


3.80 ±0.19 


9.8 ± 1.2 


-2.37 ±0.10 


-0.57 ±0.18 


0.70/6 


2008-02-11 


5.34 ± 0.32 


12 ±1 


-2.11 ±0.14 


-0.44 ± 0.24 


1.10/6 


2008-04-02 


2.94 ± 0.32 


7.1 ±0.5 


-2.44 ±0.16 




0.31/3 


2008-04-03 


4.53 ± 0.30 


11 ± 1 


-2.35 ±0.10 




0.37/6 



Table 1. Results from MAGIC-I VHE observations of Mrk 421 during the eight nights with simultaneous Swift/XKT data. 



Obs. ID 


Start 


Obs. 


f2-10keV 


«1 


E break 


"2 




^ 2 /d.o.f. 




Time (UT) 


Time 


erg/cm 2 /s 




keV 










yyyy-mm-dd hh:mm 


ks 


[•io- 12 ] 












00030352041 


2008-01-08 02:30 


2.0 


280 


2.31(0.03) 


1.20(0.09) 


2.58(0.03) 


0.242(0.003) 


1.34/163 


00030352042 


2008-01-09 04:04 


2.0 


283 


2.28(0.03) 


1.05(0.10) 


2.60(0.03) 


0.257(0.005) 


1.30/170 


00030352044 


2008-01-10 02:27 


2.3 


284 


2.32(0.02) 


1.10(0.10) 


2.57(0.02) 


0.245(0.003) 


1.59/178 


00030352053 


2008-01-16 03:21 


1.2 


345 


2.19(0.04) 


1.24(0.18) 


2.45(0.04) 


0.242(0.004) 


1.16/122 


00030352055 


2008-01-17 03:29 


0.8 


311 


2.21(0.03) 


1.97(-0.18/+0.4) 


2.75(-0.09/+0.19) 


0.243(0.003) 


1.33/101 


00030352068 


2008-02-11 03:40 


1.9 


587 


2.19(0.01) 


2.43(0.2) 


2.57(0.06) 


0.372(0.002) 


1.76/215 


00030352083 


2008-04-02 00:42 


0.9 


474 


2.09(0.02) 


2.86(0.33) 


2.51(-0.08/+0.15) 


0.260(0.003) 


1.47/130 


00030352086 


2008-04-03 21:59 


1.2 


961 


1.95(0.02) 


2.37(-0.16/+0.28) 


2.33(0.06) 


0.438(0.003) 


1.48/210 



Table 2. Mrk 421 soft X-ray fluxes and spectral parameters from the eight Swift/XKT datasets simultaneous to MAGIC-I observa- 
tions. 
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Fig. 4. VHE SED of Mrk 421 derived from the MAGIC-I observations performed in the eight time slot s with tightly simul t aneou s 
Swift/XRT data (see Section|4|. The spectra are shown after correction of the EBL absorption, following lFrance schini et al. (2008). 
For comparison, the spectrum derived from the observation that registered the highest flux (3.6 C.U. above 200 GeV) of the whole 
campaign, performed on 2008 March 30. For each night we also plot our model (see Section l4~TT > to guide the eye. 



Night 


MAGIC Obs. 


Swift/XRT Obs. 


Overlap 




Start (UT) 


End (UT) 


Time 


ZA range 


Start (UT) 


Time 


Time 


yyyy-mm-dd 


(hh.mm) 


(hh.mm) 


(ks) 


(deg) 


(hh.mm) 


(ks) 


(ks) 


2008-01-08 


01.58 


02.44 


2.7 


31-41 


02.30 


2.0 


0.8 


2008-01-09 


03.56 


06.19 


8.2 


6-20 


04.04 


2.0 


2.0 


2008-01-10 


02.23 


06.05 


3.5 


10-34 


02.27 


2.3 


1.1 


2008-01-16 


03.17 


05.11 


6.5 


6-21 


03.21 


1.2 


1.2 


2008-01-17 


03.26 


04.25 


3.4 


6-18 


03.29 


0.8 


0.8 


2008-02-11 


03:33 


03:58 


1.4 


10-18 


03.40 


1.9 


1.1 


2008-04-02 


00.45 


01.00 


0.9 


17-21 


00.42 


0.9 


0.6 


2008-04-03 


21.55 


22.20 


1.3 


15-23 


21.59 


1.2 


1.1 



Table 3. Summary of the eight tightly simultaneous observations of Mrk 421 with MAGIC-I and Swift/XRT. 
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range of each dataset. The start time and duration of the cor- 
responding Swift/XRT pointing are also reported, along with the 
actual overlapped observing time (in ks) in the last column. The 
MAGIC-I data considered for each night always cover a times- 
pan that is longer than the Swift exposures: this was necessary, 
because the typical observation time of Swift in this campaign 
(1 ks) is enough for deriving a fairly detailed X-ray spectrum of 
Mrk 42 1 , but the significantly lower count rate available in the 
y-ray domain makes this exposure time too short for deriving a 
VHE spectrum detailed enough for the modeling. Therefore the 
whole MAGIC-I exposure was used to derive the VHE spectrum 
for each night, given that the observing conditions were stable 
and no evidence for sharp evolution of the source arose from the 
VHE light curves at minute scales. For each of the eight states 
we built the MWL SED matching the MAGIC, Swift/XRT and 
optical-UV (either R-band from KVA, or UV fromSvw/f/XRT, 
or both). As an example, the SED of Mrk 421 as observed on 
2008 February 11 is plott ed in Figure |5] compared to h istori- 
cal MWL data taken from ITavecchio & Ghisellinil (120081) . It is 
worth noticing that the VHE SED is high and hard (in agreement 
with what is expected fro m other observation s of this source dur- 
ing active phases, see e.g. lAcciari et all 1201 ll) . while the X-ray 
SED is high but quite soft compared to past states where the syn- 
chrotron peak was observed at higher energies. The wide sepa- 
ration of the two peaks is discussed in detail in Section l4~T1 

4.1. SED modeling 

To reduce the degrees of freedom, we used a simp l e one- 
zone SSC model (for deta ils see ITavecchio et al 1 [19981 
Maraschi & Ta vecchiol 120031) . similar to the models com- 
monly adopted to reproduce the SEP of Mkn 421 (e.g. 
iKrawczvnski et alll2004tlFinke et all [2008). The emission zone 
is supposed to be spherical with radius R, in motion with bulk 
Lorentz factor F at an angle 9 with respect to the line of 
sight. Special relativistic effects are described by the relativistic 
Doppler factor, 6 = [F(l -Bcos8)]~ l . The energy distribution 
of the relativistic emitting electrons is described by a smoothed 
broken power law function, written for better clarity in terms of 
the adimensional Lorentz parameter y = E/m e c 2 ; the distribu- 
tion spans the [y m i n , y max ] energy range, with slopes n\ and «2 
below and above the break energy yt,, respectively. This purely 
phenomenological distribution, with n\ < 3 and n<x > 3, is able 
to reproduce the observed bumpy SED. To calculate the SSC 
emiss ion we used the full Klein-Nishina cross section dJonesL 
fl968h . 

As emphas i zed i n iBednarek & Protherod (1 19971) and 
ITavecchio et all d 1998b . constraints can be put to this simple 
model by means of simultaneous multi-wavelength observa- 
tions. Indeed, the total number of free parameters of the model 
is reduced to nine: the six parameters specifying the electron en- 
ergy distribution plus the Doppler factor, the size of the emission 
region and the magnetic field. On the other hand, from X-ray and 
VHE observations one can ideally derive seven observational 
quantities: the slopes of the synchrotron bump before and above 
the peak ai t 2 (uniquely connected to n\^), me synchrotron and 
SSC peak frequencies (v Si c) and luminosities L s c and the mini- 
mum variability timescale f var , which provides an upper limit to 
the size of the sources through the relation R < ct wm 5. It must be 
noted that as long as y m [„ <k yb <k y max the values of y^n and 
y max are not very constrained by the observation of the peaks. 
Nevertheless, the availability of the spectral shape across the in- 
strument bandpasses and of data at other wavelengths (optical 
and UV in this case) provide additional constraints with respect 



to the simple seven quantities enumerated above. Therefore, 
once all the observational quantities are known, one can fairly 
unambiguosly derive the set of parameters. In this respect, the 
cases studied here are quite favorable, because we have a fairly 
good determination of the peak frequencies (and fluxes) of both 
peaks. Indeed, although the synchrotron peak of Mrk 421 is sel- 
dom observed within the band encompassed by XRT, the joint 
optical-UV and X-ray data provide a good constraint to the po- 
sition of the synchrotron peak in all the cases. The SSC peak 
is located either within (see e.g. the SED from February 1 1 in 
Figure [3j the MAGIC-I band, or around its lower edge; in the 
latter case the pronounced curvature of the MAGIC-I spectrum 
at the lowest energies allows us to constrain the peak at energies 
not much below « 50 GeV 

Unfortunately, for the epochs used to derive the SEDs we 
do not have information on the variability timescale, f var , one of 
the key observational parameters needed to completely close the 
system and uniquely derive the parameters. In the X-ray band 
Swift/XRT observed in short (~ 1 ks) snapshots, while no con- 
vincing evidence for sub-hour variability arose from the corre- 
sponding MAGIC-I observations. Therefore we still have some 
freedom in choosing the input parameters: one can obtain dif- 
ferent sets of parameters that reproduce the spectral data equally 
well but differ in the predicted observed minimum variability 
timescale. 

We applied the model to all eight sets of data collected when 
Swift and MAGIC-I could observe the source simultaneously. 
UVOT and KVA data were also included in the SED when avail- 
able. The sets of parameters obtained from the modeling are re- 
ported in Tableland the SED data and the corresponding model 
are plotted in Figure [6] The table reports for each night the min- 
imum (y m i„), break (y/,) and maximum (y max ) Lorentz factors of 
the electron distribution, the low (n\) and high («2) energy slopes 
of the electron distribution, the magnetic field (B) and the elec- 
tron normalization (K) within the emitting region, the radius of 
the emitting region {R) and its Doppler factor (6). From these 
input parameters we derived the light crossing time (t var ); the 
contributions to the total jet luminosity from cold protons (LE . ) 
and relativistic electrons (LL) in the jet, and from magnetic field 
(Lb)', the electron (u e ) and magnetic (ug) energy densities. We 
note that in reproducing the SED we did not consider the radio 
data, since the modeled region is opaque at these frequencies: 
in this framework the radio emission originates in regions of the 
jet farther away from the black hole, beyond the core visible at 
VLBI scale, which is thought to mark the radio "photosphere". 
Accordingly, the inferred source radius is well within the up per 
limit of 0.1 pc (3 x 10 17 cm) imposed bv lCharlot etall d2006l) for 
the projected size of the SSC zone, based on VLBI observations 
of the radio core. 

Inspection of Table [4] shows that the derived Doppler factors 
are quite high, exceeding 6 — 40 in all the cases and reaching 
values as high as 80-85 in the most extreme cases. The main 
reason for these high values of 5 is the large separation be- 
tween the two peaks, the synchrotron one located below 10 17 
Hz, the SSC one around 10 25 Hz or above. As detailed in, e.g., 
ITavecchio & Ghisellinil d2008l) . a large distance between the two 
peaks implies a fairly high value of the Lorentz factor at the 
peak, yb, since yb = (vc/v s )^ 2 , and this directly implies a low B 
and a large 5 to satisfy the other constraints. 

However, we recall that because the variability timescale is 
not known, we are left with some freedom in selecting the in- 
put parameters. In the models reported in Figure |6]we assumed 
variability timescales in the range 0.5-2 hours, as typically de- 
rived for these sources (see discussion below). In general, the 
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Fig. 5. Example MWL SED of Mrk 421 as observed in tight simultaneity by Swift/UVOT, Swift/XRY and MAGIC-I on 2008 
Febr uary 1 1 (filled red circle s). Historical data (observed in different campaigns, and under less demanding time constraints) taken 
from Tavecchi o et al.l (1201 Ol) are plotted for comparison, with gray open symbols. 
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6.6 
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0.47 
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8.0 


3.2 


1.0 


2.0 


3.5 


0.050 


5900 
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70 


0.5 
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80 


0.46 
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9.2 
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17 


20 
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2.0 


4.0 


0.040 


2000 


8.5 


40 


2.0 


5.47 
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0.62 
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Table 4. Input model parameters and derived physical quantities for each of the eight simultaneous SED. 



required Doppler factor roughly scales with the observed vari - 
ability timescale as 6 oc f~£- 5 (e.g. Tav ecchio & G hisellini 2008). 
Therefore, relaxing the condition on ? var and allowing longer 
minimum variability timescales, one obtains lower 5. As an ex- 
ample we used the case for which we derive the largest 6, that of 
February 11, requiring 6 = 85. As noted above, here the deter- 
mination of the peak frequencies is very robust, because the SSC 



peak falls well within the band covered by MAGIC-I. Therefore 
this is also the best "benchmark" available to test the robust- 
ness of the derived parameters. For this purpose we modeled 
this SED assuming two sets of parameters, basically differing 
for the value of the Doppler factor, the radius of the emitting 
region and the magnetic field intensity. For 6 — 85 we have 
fvar = 0.7 h (2.5 x 10 3 s), while more than halving the Doppler 
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Fig. 6. SED of Mrk 421 with the SSC model overplotted for each of the eight simultaneous sets of MAGIC-I, Swift/XRT and 
optical-UV data obtained in the 2008 campaign. Observation dates are January 8, 9 and 10 (a), January 16 and 17 (b), February 1 1 
(c) and April 2 and 3 (d). 



factor 5 — 40 implies a fairly long variability timescale, f var = 5 
h (1.8 x 10 4 s), already longer than the characteristic variabil- 
ity timescale of Mrk 421 in the X-ray band. We can conclude 
that for the case of February 11, although the parameters cannot 
be uniquely fixed, the required Doppler factor is high, at least 
higher than 6 ~ 40. All other cases are similar. The derived 
light crossing times are within the 0.5-2 hours interval. This 
hypothesis matches well the observed typical raising/decaying 
timescales of flares of Mrk 421 and similar HBLs (PKS 2155- 
304, Mrk 501), which are characterized by doubling/halving 
times of * 1 4 s (e.g. iFossati et al.ll2008t iRavasio et afll2004t 
IZhangll2002t [fanihata et al.| l2000|h with evidence for the occur- 
rence of even faster events (e.g. Gaidos et al.l[l9 96; Cui 2004). 



However, if one relaxes this assumption on the variability 
timescales, the required Doppler factors remain high. For the 
data from February 1 1 , which allow a quite robust constraint due 
to direct observation of both the synchrotron and SSC peak, this 
implies 5 > 45. In Table|4]we also report the derived powers car- 
ried by the different jet components, namely cold protons (Lf . ), 
relativistic electrons (Lt ) and magnetic field (Lb), assuming a 
composition of one cold proton per relativistic electron. Finally, 
we computed the electron (u e ) and magnetic (ub) energy den- 
sities. The jet appears to be strongly matter-dominated, as pre- 
dicted in the standard picture of HBL sources, and i n good agree- 
ment w ith the result of the modeling performed in lAcciari et alj 
(l20Tll) . 
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5. Discussion 

During the 2008 campaign on Mrk 421 with MAGIC-I a very 
interesting dataset was gathered in VHE y-rays, complemented 
by crucial data in optical-UV and soft X-rays. For the first time 
it was possible to collect data in these bands in close simultane- 
ity during high states of the source, so that the derived spec- 
tra sampled the SED close to the synchrotron and IC peaks. In 
this situation the parameters describing the source in the frame- 
work of the standard one-zone leptonic model can be determined 
with remarkable robustness. One of the most relevant results of 
our analysis is that to reproduce the observed SED with this 
model, very high Doppler factors are required. There is some 
freedom in choosing the parameters, mainly because of the not 
known variability timescale at those epochs. In the models sum- 
marized in Table |4] and reported in Figure |6] we assumed vari- 
ability timescales in the range of 0.5-2 hours. 

Indeed, these high values of inferred 5 are not rare: 
very high Doppler factors, sometimes higher than 5 ~50, 
for Mrk 421 and other well-observed HBLs were ob- 
tained in the past, leading to the so called "<5-crisis" (e.g . 
Krawczvnski et alll2002tlAcciari et all 1201 \\ [Konopelko et all 
2003; Georganopoulos & Kazanasl 120031: iHenri & Sauga, 
2006; iGiebels et all 12001 iFinke et all l2008h . Analogously. 
the recent exc e ptiona l VHE flare of PKS 2155-304 
( Aha ronian et all 120071) seems to require e xtreme Doppler 
factor s i n the framework of one -z one m o dels dBegelman et aD 
20081 iGhisellini & Tavecchid 120081 IFinke et all 120081: 
Kusunose & Takaharal 20061) . These high values of 5 (im- 
plying a similarly high value of the bulk Lorentz factors) 
contrast with the very low jet velocities infe rred at VLB! 
scales in a large part of Te V BL Lacs (e.g. iGiro letti et al. 
120041 iPiner & Edwardsl 120041) . including Mrk 421, and with 
the value of T required fr om the unification of BL L a cs and 
FRI radiogalaxies ( e.g. iGeorganopoulo s & Kazanasl 120031: 
lHenri&Sau"gl.l2006h . 

In addition to the extreme Doppler factor, we can identify 
two other problems afflicting the standard model, namely the 
huge difference between the magnetic and particle energy den- 
sities and the extremely long cooling times of the emitting rela- 
tivistic electrons. 

Table [4] shows that in all cases the inferred electron energy 
density substantially exceeds the corresponding magnetic energy 
density, by up to two orders of magnitude and even more. This 
resul t is also generally f ound from SSC fits of HBL SEDs (see 
e.g. lAcciari et all 1201 ll for a recent example), while for other 
classes of blaza rs, in particular for FSR Qs, equipartition is usu- 
ally found (e.g. IGhisellini et all 12010). This evidence disagrees 
with the general expectations of the diffusive shock acceleration 
models, in which a substantial equipartition between particles 
and magnetic field is expected. 

Concerning the second problem, namely the long cooling 
times cales, we point out that, follow in g for instance the formu- 
lae in lBednarek & Protheroel (1 19971) or ITavecchio et al.l (Tl998). 
fairly long cooling times f coo /, on the order 10 6 s in the observer's 
frame, can be computed from the model parameters in Table 2] 
Therefore, the declining part of flares cannot be attributed to the 
cooling of the emitting electrons. A possibility is that adiabatic 
expansion, which allows quenching of the flux within scales of 
R/c, has to be invoked as one of the viable processes that may 
explain the observed descent of TeV and X-ray fluxes on hour 
scales. But it must be noted that this explanation has the signifi- 
cant drawback of implying a very energetically inefficient jet. A 
problem related to the long cooling timescale is that one cannot 



interpret the observed break in the underlying electron energy 
distributio n as the separation betw een fast and slow cooling par- 
ticles (e.g. ITavecchio et all U998). so it has to be, for instance, 
assumed to be intrinsic to the injected population. These prob- 
lems could therefore hint at the unsuitability of the one-zone 
model for this source. 

A solution of these problems faced by the standard one-zone 
SSC scenario, extensively discussed in literature, is based on 
the possible existence of multiple active emitting regions (e.g . 
Blazeiows ki et all 120051 : IGeorganopoulos ^ Kazanasl [2003). 
based on the possibility that the flow is characterized by por- 
tions moving at different speeds. If these regions emit, in each of 
them the electrons can scatter not only the locally-produced syn- 
chrotron photons, but also the soft photons produced in the other 
region. Moreover, the energy density of these "external" photons 
is amplified in the rest frame of the emitting region through the 
relative speed between the two portions of the flow. The result 
is that the inverse Compton emission of each region is ampli- 
fied with respect to the SSC emission. As a consequence, the 
Doppler factor required to reproduce the SED is lower than that 
of the one-zone model. I n particular, in the "spine-layer" model 
of Ghisell ini et alJ d2005l) . it is assumed that the jet has an inner 
faster core (the spine) that is surrounded by a slower layer. At 
a narrow angle of view, which is characteristic for blazars, the 
emission is dominated by the faster spine whose IC emission is 
a mixture of SSC and "external" Compton components. 

This model would also more easily accomodate the short 
variability time scales observed in Mrk 421, which are hardly 
explained within the one-zone model due to the long electron 
cooling times, as discussed in Section |4~T1 Indeed, in general a 
10 times more intense B can be adopted when modeling a given 
SED (see e.g. IGhisellini et al]|2005l) : because the synchrotron 
cooling time scales as f JV „ c oc B 2 , this could lead to cooling 
times on the order of the required variability time scale. 

An alternative scenario that is possibly able to solve these 
problems i s the "minijets" model advocated by iGiannios et al.l 
(l2009Ll2010h . In this framework, the emission is thought to occur 
in very fast small portions of plasma resulting from the rapid 
reconnection of magnetic field lines inside the main jet flow. If 
the magnetization (ratio of magnetic over kinetic jet power) is 
high enough, the Lorentz factor of these blobs in the rest frame 
of the jet can be as high as F = 50. Moreover, since the emitting 
plasma is the residual of the annihilation of magnetic field, one 
naturally expects a low magnetic energy density and thus a high 
particle over magnetic energy density ratio. 

A modeling of the SED with the more complex (and less 
constrained) models mentioned above is beyond the scope of this 
paper and left to future work. 
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